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Energy levels of germanium, tin and lead together with their single, double and triple ionized 
positive ions have been calculated using the y^~*^ approximation suggested in the previous work 
(Dzuba, physics/0501032 I (M = 4 - number of valence electrons). Initial Hartree-Fock calculations 
are done for the quadruply ionized ions with all valence electrons removed. The core-valence cor- 
relations are included beyond the second-order of the many-body perturbation theory. Interaction 
between valence electrons is treated by means of the configuration interaction technique. It is demon- 
strated that accurate treatment of the core-valence correlations lead to systematic improvement of 
the accuracy of calculations for all ions and neutral atoms. 
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I. INTRODUCTION 

This work further develops the V^~^^ approximation 
suggested in Ref. Ij. It also presents the details of the 
calculations of the energy levels of Ge II, Sn II and Pb II 
which were needed to study their dependence on the 
fine structure constant a (a = e^/fic). Some lines of 
Ge II, Sn II and Pb II have been observed in quasar ab- 
sorption spectra and the information on the dependence 
of corresponding frequencies on a is needed to study 
possible variation of the fine structure constant at early 
epoch. 

In the vicinity of the physical value of a the frequency 
of an atomic transition can be presented in a form 

U} = uJo + qx, (1) 

where x = {a/aof' — 1 and and wq are the present- 
day laboratory values of the fine structure constant and 
transition frequency. 

The values of the g-coefficients can only be found from 
atomic calculations by, e.g., varying the value of a in 
computer codes based on relativistic equations. In many 
cases calculated values of the g-coefiicients are more sta- 
ble than the energies. This is because they are not sensi- 
tive to incompleteness of the basis set with respect to the 
principal quantum number n. Indeed, relativistic correc- 
tions are proportional to Q {v is the eff'ective prin- 
cipal quantum number) while energies are proportional 
to If we include more states of high v this would 

have greater effect on the energies than on relativistic 
corrections presented by g-coefficients. 

However, in the case of strong configuration mixing 
and level pseudo-crossing calculation of g-coefficients may 
become very unstable . In the vicinity of level pseudo- 
crossing the values of g-coefficients change very rapidly 
with a and small error in determining the position of the 
level crossing may lead to large error in the values of q. 
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Level pseudo-crossing always means strong configura- 
tion mixing between the states. However, strong config- 
uration mixing may also take place without level pseudo- 
crossing. This can also cause instability in calculated 
values of ^-coefficients. Indeed, relativistic correction to 
the energy of a single electron state \njlm) strongly de- 
pends on the total momentum j of this state (see, e.g. 
formula (7) in Ref. [sf). Therefore configurations com- 
posed from states of different j may have very different 
values of q and small error in the the configuration mix- 
ing coefficients would lead to large error in the resulting 
q value for the mixed state j^. 

Strong configuration mixing and level pseudo-crossing 
take place for Ge II, Sn II and Pb II ions ^ as well as for 
many other atoms and ions This means that calcu- 
lations need to be done to very high accuracy to ensure 
stable values of the g-coeflacients. The criterion is that 
deviation of the calculated energies from the experimen- 
tal values must be much smaller than the experimental 
energy interval between mixed states. 

There are many other areas of research where accurate 
atomic calculations are needed. These include parity and 
time invariance violation in atoms (see, e.g. [Sl), atomic 
clocks 7j, interaction of positrons with atoms [8j, etc. 

A way to do accurate calculations for atoms with sev- 
eral s and/or p valence electrons has been suggested in 
Ref. 0. It is called "the y^^*^ approximation", where 
V is the Hartree-Fock potential created by — M elec- 
trons of the closed shell ion, A^ is total number of elec- 
trons in neutral atom and M is the number of valence 
electrons. Initial Hartree-Fock calculations are done for 
a closed-shell positive ion with all valence electrons re- 
moved. It has been demonstrated in Ref. j3| that the 
Hartree-Fock potential of the closed-shell positive ion is 
often a good starting approximation for a neutral atom. 
This is the case when valence electrons are localized on 
distances larger than the size of the core. Then they 
can affect only energies of core states but not their wave 
functions. Since the potential created by core electrons 
depends on the electron charge density and does not de- 
pend on electron energies it doesn't matter which core 



2 



states are used to calculate the potential - states of the 
neutral atom or states of the closed-shell positive ion. 

The effective Hamiltonian for valence electrons is con- 
structed using the configuration interaction (CI) tech- 
nique. Core- valence correlations are included by adding 
the electron correlation operator S to the CI Hamilto- 
nian. Many-body perturbation theory (MBPT) is used to 
calculate S. The main advantage of the v^~^'^ approxi- 
mation is that MBPT is relatively simple (no subtraction 
diagrams) and the S operator can be calculated beyond 
the second-order of the MBPT. It has been demonstrated 
in Ref. 1] that inclusion of the higher-order core valence 
correlations lead to further significant improvement of 
the accuracy of calculations. 

In the previous work 1] the v^^^'^ approximation was 
used for Kr and Ba while higher-order core-valence cor- 
relations were included for Ba and Ba"*" only. In the 
present work we study twelve complicated many-electron 
systems including germanium, tin, lead and their positive 
ions. We demonstrate that using the V^~'^ approxima- 
tion {M ~ 4 for the case of Ge, Sn and Pb) and accurate 
treatment of the core-valence correlations lead to high 
accuracy of calculations for all twelve systems. This in- 
dicates that the v^~^^ approximation is a good approx- 
imation for a wide range of atoms and ions. 

II. CALCULATIONS 

The effective Hamiltonian for valence electrons in the 
yN-M appj-Qximation has the form 

A I M 

H-^ =Y,hu+Y,h2,„ (2) 

i=l i^j 

hi{ri) is the one-electron part of the Hamiltonian 

/ii =ca-p+(/3-l)TOc2- — + y^-4 + Ei. (3) 

r 

Si is the correlation potential operator which is exactly 
the same in the V'^~^^ approximation as for the single- 
valence electron atoms (see, e.g. 0). It can be calculated 
in the second-order of the MBPT. Selected chains of the 
higher-order diagrams can be included into Ei in all or- 
ders using technique developed for single- valence electron 
atoms (see, e.g. [13 )• 

h2 is the two-electron part of the Hamiltonian 

/l2 = l r+S2(ri,r2), (4) 

In -1-21 

is the two-electron part of core-valence correlations. 
It represents screening of Coulomb interaction between 
valence electrons by core electrons. We calculate S2 in 
the second order of MBPT. Inclusion of the higher-order 
correlations into E2 will be a subject of further study. 
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FIG. 1: Electron density of the s,p,d and / electrons of Pb I 
and Pb V as explained in the text. 



However, the calculations show that in most cases accu- 
rate treatment of Si is more important than for S2. The 
details of the calculation of Si and S2 can be found else- 
where (9l il0llTll [l^. Note however that in contrast to the 
previous works ^3 0| we have no so called subtraction 
diagrams. 

Number of electrons M is the only parameter in the ef- 
fective Hamiltonian (|2J) which changes when we move be- 
tween different ions of the same atom. The terms V^~'^, 
Si and S2 remain exactly the same. 

The form of the effective Hamiltonian is also the same 
for all ions if some other potential V is used to gener- 
ate the core states. However, the S operator would have 
terms proportional to V^~'^ — V (subtraction diagrams 
[ni). In the T/JV-M approximation V = l/""" and sub- 
traction diagrams disappear. The MBPT becomes rel- 
atively simple which makes it easier to include higher- 
order core- valence correlations. 



A. Electron shell structure of lead. 

To understand how the v^~^'^ approximation works 
it is very instructive to look at electron shells of a many- 
electron atom. We chose lead because it is the heaviest 
of the considered atoms. It probably has the richest pos- 
sible electron shell structure. Neutral lead has eighty 
two electrons occupying six shells. Angular momentum 
I ranges from (s-electrons) to 3 (/-electrons). Figs. ^ 
and [3 present electron densities of Pb I (solid line) and 
Pb V (dotted line) separately for s, p, d and / electrons. 
The density is the sum over principal quantum number 
n, total momentum j and its projection m while angular 
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momentum I is fixed: 



njm 



The values of p{r)i in the maximum are very different for 
different I. Therefore, we present normaUzed functions 
pii^)i/ Pmax to be able to fit all graphs into one diagram. 

Electron shell structure can be clearly seen on Fig. ^ 
Each density has a local peak aX n — I = 1,2, etc. The 
position of the peak depends mostly on n and is about 
the same for all I. This means that all electrons with 
the same n are localized at about the same distances 
regardless of their angular momentum I, thus making a 
shell. 

The difference between Pb I and Pb V cannot be seen 
on Fig. n Fig. |21 presents details of the right bottom 
corner of the Fig. ^ Dotted lines which correspond to 
electron densities of the Pb V ion show no peak at n — 
6 because of absence of the 6s and 6p electrons. The 
removal of four valence electrons has some effect on the 
density of d-electrons at about the same distances where 
the 6s and 6p electrons are localized and practically no 
effect on the densities of all electrons on shorter distances. 
This is because valence electrons are localized on large 
distances and they can only create constant potential in 
the core which can change the energies of the core states 
but cannot change their wave functions. 

One can see from Fig. |5] that there is an overlap be- 
tween the wave functions of valence electrons of Pb I (6s 
and 6p electrons) and the wave function of the core out- 
ermost state 5d. We have presented for comparison on 
Fig. O the electron densities of Ba I and Ba III on large 
distances. It is easy to see that the overlap between core 
and valence electrons in barium is much smaller than the 
overlap between core and valence electrons in lead. As 
a consequence, outermost core state of barium (5p) is 
much less affected by removal of two 6s electrons than 
compared to the effect of removal of two 6s and two 
6p electrons on the 5d state of lead. This means that 
the V^^"^ approximation for Ba should work much bet- 
ter than the V^^^ approximation for Pb. The situation 
is exactly the same as for the V^~^ approximation for 
atoms with one external electron. It is very well known 
that the V^~^ approximation works extremely well for 
alkali atoms and not so well for atoms like Ga, In, Tl, etc. 
The reason is the same in both cases. Valence electrons 
must not overlap with the core for the V^~^ to be good 
starting approximation regardless of whether M = 1 or 
M > 1. 

Similar to the fact that the V^~^ approximation is a 
good approximation for thallium, although not as good 
as for alkali atoms, the approximation is a good 

approximation for Pb, Sn and Ge, although not as good 
as for Ba. 

Below we present specifics of calculations for gcrma- 
nimn, tin and lead. 
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FIG. 2: Details of electron densities of Pb I and Pb V at large 
distances. 




5 10 
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FIG. 3: Electron densities of Ba I and Ba III at large dis- 
tances. 



B. Calculations for germanium 

Germanium is the lightest of three atoms (Z=32) and 
the easiest from computational point of view. Its ground 
state configuration is ls^2s^2p^3s^3p'^3d^°4:s'^ip'^ . The 
core- valence correlations are relatively small due to small 
number of electrons in the core. 

We calculate Ei and S2 for the effective Hamiltonian 
(0) in the second order of the MBPT. Inclusion of Si 
brings single-electron energies of Ge IV to agreement with 
the experiment on the level of 0.1%. No higher-order 
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core-valence correlations need to be included. 

In fact, inclusion of the higher-order correlations using 
technique developed in Ref. doesn't lead to better 
results for germanium. This is because the technique 
was developed for heavy atoms in which higher order 
correlations are dominated by screening of the Coulomb 
interaction between core and valence electrons by other 
core electrons. In light atoms like germanium this effect 
does not dominate due to small number of electrons in 
the core. Therefore, inclusion of screening, while other 
higher-order effects are not included, does not improve 
the accuracy. 

The results of calculations are presented in Table 
The ground-state energies are given as energies to re- 
move all valence electrons from an atom or ion (in atomic 
units). Corresponding experimental energies are sums 
of the ionization potentials of all relevant ions. For the 
convenience of comparison with Moore's tables '^6 
present energies of excited states relative to the ground 
state in cm~^. Column marked CI presents the results of 
the standard configuration interaction method without 
S. Column E^^^ presents the results of calculations with 
the effective Hamiltonian ^ in which E is calculated in 
the second order of MBPT. 

The results presented in Table show that inclusion of 
the core-valence correlations leads to systematic signifi- 
cant improvement of the accuracy of calculations for all 
states of all ions and for neutral germanium. 

C. Calculations for tin. 

Tin atom {Z = 50) is very similar to the germanium 
atom. Its ground state configuration is . . . 5s^5p^. How- 
ever, correlations and relativistic corrections are larger. 
It has some implication on the calculation scheme. It 
turns out that inclusion of the higher-order core-valence 
correlations does lead to significant improvement of the 
results for all tin ions and for the neutral atom. We in- 
clude screening of Coulomb interaction and hole-particle 
interaction in all orders of the MBPT in the calculation 
of El. It is done exactly the same way as in our calcu- 
lations for single- valence-electron atoms (see, e.g. [Tojl. 
The E2 operator is still calculated in the second order of 
the MBPT. 

The results are presented in Table ^ There is one 
more column in the table compared to Table P It is 
marked E(°°) and presents the results of calculations with 
all-order Ei. Again, it easy to see that moving from 
less sophisticated to more sophisticated approximations 
(with no E; with E^^'; with E(°°)) leads to systematic 
significant improvement of the accuracy of the results. 

D. Calculations for lead 

The case of lead {Z = 82) is the most difficult of the 
calculations. Correlations are strong and relativistic ef- 



TABLE I: Ground state removal energies (RE, a.u.) and ex- 
citation energies (cm~^) of low states of Ge IV to Ge I. 
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fects are large too. Strong L — S interaction leads to 
intersection of the fine-structure niultiplets. Also, states 
of the same total momentum J are strongly mixed re- 
gardless of the values of L and S assigned to them. The 
breaking of the L — S scheme can be easily seen e.g. 
by comparing experimental values of the Lande g-factors 
with the non-relativistic values. 

We have done one more step for lead to further im- 
prove the accuracy of calculations as compared to the 
scheme used for tin. We have introduced the scaling fac- 
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TABLE II: Ground state removal energies (RE, a.u.) and 
excitation energies (cm~^) of low states of Sn IV to Sn I. 

State CI S'^^ S'°°' Exp. [14] 

Sn IV 



tors before Si to fit the energies of Pb IV. These energies 
are found by solving Hartree-Fock-hke equations for the 
states of external electron of Pb IV in the ^^-^ potential 
of the atomic core 

[Ho + Si - e,Mn = 0. (5) 

Here is the Hartree-Fock Hamiltonian. Ei is the 
all-order correlation potential operator similar to what is 
used for tin. Inclusion of Ei takes into account the effect 
of the core- valence correlations on both the energies (e„) 
and the wave functions (V'n) of the valence states pro- 
ducing the so-called Brueckner orbitals. The difference 
between Brueckner and experimental energies of the 4s, 
4p and 4c? states of Pb IV are on the level of 0.2 - 0.4% 
(for removal energies). To further improve the energies 
we replace Ei by /Ei with rescaling factor / chosen to 
fit the energies exactly. Then the same rescaled operator 
/El is used for the Pb III and Pb II ions and for the 
Pb I. It turns out that only small rescaling is needed. 
Maximum deviation of the rescaling factor from unity 
is 10%: /(4s) = 0.935, /(4pi/2) - 1.084, f{Ap^,2) = 
1.1, /(4d3/2) = 1.07, /(4d5/2) = 1.07. 

The results of the calculations are presented in Ta- 
ble Again, inclusion of core- valence correlations lead 
to significant improvement of the accuracy of the results 
in all cases. However, comparison between different ways 
of treating core-valence correlations reveal a more com- 
plicated picture compared to what we have for tin. When 
we move from the second-order correlation operator E'^' 
to the all-order E(°°) and then to the rescaled /S(°°) 
the improvement in accuracy is apparent for the removal 
energies. It is again systematic and significant, bring- 
ing results for all states of all ions and neutral lead to 
better agreement with experiment. This is not always 
the case for the energy intervals. When a more accurate 
treatment of core-valence correlation is introduced two 
energy levels way move cowards experimental values at 
slightly different rate so that the interval between them 
does not improve. In Table UTTI we present removal ener- 
gies only for the ground states of Pb IV, Pb III, Pb II 
and Pb I. Energies of excited states are given with re- 
spect to the ground state. It is easy to see that energy 
intervals between ground and excited states calculated 
with second-order E are often in better agreement with 
experiment than the results with the all-order El. In gen- 
eral, the results are not as good as for tin. The reason 
for this is larger overlap between valence and core states. 
Relativistic effects cause stronger binding of the 6s and 
6p electrons of Pb compared to binding of the 5s and 5p 
electrons of Sn. This means that overlap between valence 
and core states is also larger for lead than for tin leading 
to larger effect of removal of valence electrons on atomic 
core. 

It is instructive to compare our results with the re- 
sults of recent calculations by Safronova et al (see 
Table Unil. Energy levels of Pb II were calculated by 
Safronova et al with the use of the coupled-cluster (CC) 
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approach and the thhd-order MBPT. The Pb II ion was 
treated as an ion with one external electron above closed 
shells. Therefore only energies of states in which the 6s 
subshell remained closed were calculated. The agreement 
with experiment for these states is slightly better than 
for our results with The reason for this is better 

treatment of the interaction between core and valence 
electrons. The 6s electrons were included in the initial 
Hartree-Fock procedure. Also, interaction between the 
6p electron and the core is included in the CC approach 
in all-orders of the MBPT. 

This doesn't mean that the V^~'^ approximation is 
not good for lead. First, as can be seen from Table ITTll 
inclusion of core- valence correlation does lead to system- 
atic significant improvement of the accuracy and final 
results are very close to the experiment. Second, the 
fact that inclusion of the higher order core-valence cor- 
relations doesn't always lead to improvement of energy 
intervals doesn't mean that the V^~'^ approximation is 
not good. It rather means that not all dominating higher- 
order diagrams are included into The situation 
is very similar to what takes place for single-valence- 
electron atoms. The technique developed by us for alkali 
atoms ^3] doesn't work very well for atoms like thallium 
where interaction between valence electron and the core 
is important. Here CC-I-MBPT approach gives better re- 
sults which may mean that the combination of the 
CC approach with the CI method is a better option for 
atoms like lead. This approach was recently considered 
by Kozlov 16] and Johnson However, no calcula- 
tions for real atoms have been done so far. 



III. CONCLUSION 



It has been demonstrated that the 1/^-* approxima- 
tion works very well for the four-valence-electrons atoms 
like germanium, tin and lead as well as for their single, 
double and triple ionized ions. The use of the y^^** ap- 
proximation makes it easy to include core- valence correla- 
tions beyond the second order of the MBPT. Inclusion of 
the core- valence correlations leads to significant improve- 
ment of the results in all cases. In general, the ]/^~*^ 
approximation (M is the number of valence electrons) is 
a good approximation if the overlap between core and va- 
lence states is small. The best case is the alkaline-earth 
atoms where the V^"'^ approximation must produce ex- 
cellent results. In contrast, the y^^*^ approximation is 
not applicable at all to atoms with open d or/and / shells 
unless uppermost core s and p states are also treated as 
valence states. It should work more or less well for most 
of the atoms/ions with s and/or p valence electrons. In 
cases of relatively large overlap between core and valence 
states good results can still be achieved if accurate treat- 
ment of the interaction between core and valence elec- 
trons is included perturbatively into the calculation of 
the core- valence correlations. 
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TABLE III: Ground state removal energies (RE, a.u.) and excitation energies (cm ^) of low states of Pb IV to Pb I. 

State CI /S<°°) Ref. [15] Exp. [14] 

Pb IV 



5d"'6s 


'S'l/a RE 


-1.48374 


-1.57689 


-1.56035 


-1.55529 




-1.55531 




'A/2 


72857 


78055 


78239 


76144 




76158 




JP3/2 


92301 


99817 


99388 


97276 




97219 




O3/2 


173446 


188501 


185992 


184570 




184559 




'O5/2 


175485 


190789 


188254 


186848 




186817 






Pb III 












^So RE 


-2.58923 


-2.76503 


-2.73356 


-2.72421 




-2.72853 


6s6p 




52866 


62881 


62947 


61045 




60397 




'Pi 


57184 


66767 


66751 


64851 




64391 




'P2 


70223 


82032 


81477 


79577 




78985 




'Pi 


91945 


96556 


95876 


94071 




95340 




'Po 


135286 


145385 


145400 


141555 




142551 


6s7s 




137664 


153445 


150863 


150038 




150084 


6s6d 




138279 


156137 


154498 


152079 




151885 


6s7s 




142139 


156815 

Pb II 


154219 


153407 




153783 




'A°/2 RE 


-3.11363 


-3.31759 


-3.27430 


-3.26897 




-3.28141 




2 po 
-^3/2 


12390 


14447 


13858 


13896 


14137 


14081 




'Pl/2 


50298 


59934 


59934 


58052 




57911 




'A/2 


57209 


68501 


67633 


66221 




66124 




JP5/2 


61484 


75957 


74856 


73749 




73905 




Sl/2 


55451 


60525 


58170 


59203 


58967 


59448 






66823 


71130 


69314 


69256 


70229 


68964 




D3/2 


63732 


70711 


68916 


69001 


69686 


69740 




2 po 
-fl/2 


69961 


75342 


73140 


73878 


74256 


74459 




2 po 
-^3/2 


72572 


78180 


75935 


76666 


77069 


77272 




O3/2 


77272 


85538 


84523 


83196 




83083 




'i?5/2 


81630 


91291 


89614 


88800 




88972 






Pb I 












'Po RE 


-3.36433 


-3.58255 


-3.53174 


-3.52974 




-3.55398 




'Pi 


6388 


7736 


7305 


7353 




7819 




'P2 


9199 


10795 


10277 


10423 




10650 






18578 


21793 


20780 


20979 




21458 




'So 


26998 


30355 


29185 


29412 




29467 


6s^6p7p 


'Po 


33413 


35239 


33679 


34517 




34960 




'Pi 


33871 


35610 


34056 


34887 




35287 


6s^6p7p 


'Pi 


40029 


42987 


41405 


42061 




42919 




'Po 


41612 


44441 


42882 


43525 




44401 


6s^6p7p 


'Ol 


41740 


44714 


43129 


43773 




44675 




'732 


41886 


44868 


43281 


43958 




44809 



